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Procede et disposi-bif d'imagerle utilisan-b des ondes de ex- 
saxllement ^ 

La presente invention est relative aux procedes et 
5 dispositifs d'imagerie utilisant des ondes de cisaillement . 

Plus particulierement, 1' invention concerne un pro- 
cede d'imagerie utilisant des ondes de cisaillement pour 
observer un milieu viscoelastique diffusant qui contient 
des particules ref lechissant les ondes ultrasonores de com- 
10 pression, procede comprenant: : 

(a) una etape d^ excitation au cours de laquelle on 
genere une onde 61astique de cisaillement dans le milieu 
viscoelastique, 

(b) une 6tape d' observation au cours de laquelle on 
15 observe la propagation de I'onde de cisaillement simultane- 

ment en une multitude de points d'un champ d' observation 
dans le milieu viscoelastique, cette etape d' observation 
comprenant les sous-etapes suivantes : 

(bl) faire emettre dans le milieu viscoelasti- 

20 que, par un reseau de transducteurs commandes ind6pendam- 
ment les uns des autres, une succession de tirs d' ondes ul- 
trasonores de compression non ,f.ac^J-±s^es- t- une ca Jerice d*' au 
moins 500 tirs par seconde, 

(b2) faire detecter et enregistrer en temps reel 

25 des signaux acoustiques regus du milieu viscoelastique 2, 
comprenant les 6chos generes par les ondes ultrasonores de 
compression non focalis^es en interagissant avec les parti- 
cules ref lechissantes dudit milieu viscoelastique, 

(c) et au moins une §tape de traitement au cours de 
30 laquelle : 

(cl) on traite les signaux acoustiques succes- 

s±fs — ^e^^u-s — ^ — milieu — vriscoelastxque au cours de la sous- 

etape (b2) pour determiner des images de propagation suc- 
cessives de I'onde de cisaillement. 
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(c2) et pour determiner au moins un parametre de 
Tnouvement du milieu viscoelastique en differents points du 
champ d^ observation . 

On obtient ainsi un "film" illustrant clairement la 
5 propagation de I'onde de cisaillement dans le milieu vis- 
coelastique^ qui peut permettre une analyse qualitative 
et/ou quantitative pour reperer des zones de durete diffe- 
rente du reste du milieu viscoelastique ou des zones ayant 
un temps de relaxation different du reste du milieu viscoe- 
10 lastique. 

Le docxament WO-A-00/55 616 d^crit un exemple d^un 
tel proced6r dans lequelon g6nere I'onde de cisaillement 
^ la surface du milieu visco61astique . Ce proced6 donne 
toute satisfaction notamment pour imager des zones situees 
15 relativement pres de la surface du milieu viscoelastique - 
Mais ce proced6 connu ne permet pas d" observer certaines 
zones du milieu viscoelastique, notamment : 

- des zones suffisamraent profondes pour ne pas 
pouvoir etre atteintes par des ondes de cisaillement gene- 

20 rees en surface (les ondes de cisaillement s'attenuent ra- 
pidement ) , 

- et des zones d' ombre masquees par des obstacles 
(notamment des parties du squelette d'un patient ou des zo- 
nes liquides telles que des kystes liquides) qui g&nent la 

25 propagation des ondes de cisaillement. 

De plus, si le champ d' observation est partielle- 
ment dans une zone d' ombre ^ il peut §tre necessaire de de- 
placer le dispositif generateur des ondes de cisaillement 
en cours d' observation, ce qui est fastidieux pour I'utili- 

30 sateur* 

Enfin, le dispositif generateur des ondes de ci- 
saillement est relativement lourd et complique le disposi- 
tif- 

La presente invention a notamment pour but de pal- 
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lier ces inconv^nients • 

A cet effet, selon 1' invention, un proc^de du genre 
en question est caract6ris6 en ce qu'au cours de l''6tape 
d' excitation (a), on fait gen6rer • 1'' onde elastique de ci~ 
5 sailleiaent en faisant emettre au moins une onde ultrasonore 
focalisee dans le milieu viscoelastique par ledit reseau de 
transducteurs, la focalisation et la chronologie de ladite 
onde ultrasonore focalisee, ainsi que la chronologie desdi- 
tes ondes ultrasonores non focalisees, etant adaptees pour 
10 qu'au moins certaines desdites ondes ultrasonores non foca- 
lisees parviennent dans le champ d'' observation lors de la 
propagation de I'onde de cisaillement dans ce champ 
d' observation . 

Ainsi, le mSme reseau de transducteurs permet ^ la 
15 fois d'engendrer I'onde elastique de cisaillement de fagon 
choisie dans le champ d' observation, et d'' observer ensuite 
cette propagation, grSce au fait que le dispositif 
d' imagerie est adapte pour generer soit des ondes ultraso- 
nores focalisees permettant de generer I'onde elastique de 
20 cisaillement, soit des ondes ultrasonores non focalisees 
permettant d' observer la propagation de I'onde de cisaille- 
ment, et gr§ce au choix judicieux : - . - 

- de la chronologie de ces differentes emissions, 
et du ou des points de focalisation de I'onde 
25 ultrasonore focalisee. 

Le precede d' imagerie selon 1' invention est done 
ais6 ^ mettre en oeuvre pour un utilisateur, avec un dispo- 
sitif relativement simple et leger. L' invention est done 
d'un tres faible cout de revient en comparaison de techni- 
30 ques concurrentes telles que I'^IRM, et elle permet le cas 
echeant de constituer des systemes d' imagerie ambulatoires 
utilisables ta nt — en — ima goric p^ upe xa foire qu'en imagerie 
postoperatoire, voire en imagerie per-operatoire . 

A titre d'exemple, dans des applications medicales. 
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le proc6d6 selon 1' invention peut permettre de reperer ef- 
ficacement des zones canc6reuses dans les tissus d'un pa- 
tient- La propagation des ondes de cisaillement s'y deroule 
en effet tres dif f erenmient des zones voisines. Ce reperage 
5 s'effectue beaucoup plus facilement que par une observation 
classique par simple echographie ultrasonore^r puisque la 
propagation des ondes de cisaillement est fonction du mo- 
dule de cisaillement du milieu, lui-meme tres variable en- 
tre une zone de tissus sains et une zone de tissus canqe- 

10 reux : le module de cisaillement varie typiquement dans un 
rapport de 1 a 30 entre une zone saine et une zone cance- 
reuse, alors que le module de compression, qui regit la 
propagation des ondes acoustiques de compression utilis6es 
dans 1 ' echographie ultrasonore, varie seulement de I'ordre 

15 de 5% entre un tissu sain et un tissu canc6reux. 

De m§me, on peut ainsi rep§rer des zones necrosees 
au sein d'un tissu^ par exemple des zones tumorales ayant 
subi un traitement par hyperthermie aux ultrasons, notam- 
ment pour evaluer I'efficacite de ce traitement par hyper- 

20 thermie, 

Une autre application possible de 1' invention 
concerne 1' Evaluation quantitative du degre de fibrose du 
foie, qui est un parametre important dans les pathologies 
du foie, notamment 1' hepatite C. 

25 On notera que 1' invention permet de generer I'onde 

de cisaillement, et d' observer sa propagation, y compris a 
travers une zone liquide ou une barriere osseuse masquant 
totalement ou partiellement le champ d' observation (crane, 
cage thoracique, etc), puisqu' il est possible de focaliser 

30 des ondes ultrasonores au travers de telles barrieres (voir 
notamment le document WO-A-02/32316 ou la demande de brevet 
franq:ais n^ 02 10682 du 28 aout 2002) . 

Dans des modes de realisation preferes du precede 
selon 1' invention, on peut eventuellement avoir recours en 
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outre a I'une et/ou a 1' autre des dispositions suivantes : 

au cours de la sous-etape (b2) pour determiner 
ledit parametre de mouvement, on compare plusieurs images 
successives de propagation (par exemple par correlation^ 
5 Doppler, etc.) avec una m§me image de reference du milieu 
viscoelastique, cette image de ref§rence 6tait determinee 
en tirant au moins une onde ultrasonore de compression non 
focalisee dans ledit milieu viscoelastique puis en d^tec- 
tant et en enregistrant des echos g6ner6s par ladite onde 

10 ultrasonore de compression non focalisee lorsqu'elle inte- 
ragit avec les particules r6f 16chissantes du milieu viscoe- 
lastique (on ameliore ainsi la precision de la mesure du 
parametre de mouvement (par exemple le deplacement) du mi- 
lieu viscoelastique, notamment pour de faibles amplitudes 

15 de mouvement (typiquement moins de 30 |jm avec la technique 
d' excitation par ondes acoustiques utilis6e ici) ; 

- I'etape (a) est precedee d'une §tape observa- 
tion initiale (aO) au cours de laquelle on tire au moins 
une onde ultrasonore de compression non focalisee puis on 

20 detecte et on enregistre des echos gen6res par ladite onde 
ultrasonore de compression non focalisee en interagissant 
avec les particules ref lechissantes du m.i lieu vi-GCOclatici- 
que, ces echos correspondant (directement ou indirectement) 
a une image initiale du milieu viscoelastique, et au cours 

25 de la sous-etape (b2) ladite image initiale constitue la- 
dite image de reference pour traiter au moins certaines des 
images successives de deplacement ; 

au cours de l''etape d'' observation initiale (aO) , 
on tire successivement plusieurs ondes ultrasonores de com- 

30 pression non focalisees puis on detecte et on enregistre 
des echos generes par chaque onde ultrasonore de compres- 
sion non focall _sep pn-^r^er-atrirs^sHrrb" avec les part i nnl p*.c;s t- e - 
f lechissantes du milieu viscoelastique, ces echos corres- 
pondant (directement ou indirectement) a plusieurs images 
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preliminaires- successives du milieu viscoelastique, et on 
determine ladite image initiale du milieu viscoelastique en 
combinant lesdites images preliminaires successives ; 

- ledit parametre de mouvement est un dfeplacement 
5 du milieu viscoelastique ; 

l''onde ultrasonore focalis6e emise au cours de 
I'etape d'' excitation (a) pr^sente une. frequence f comprise 
entre 0,5 et 15 MHz et est emise pendant une duree k/f en 
secondes, ou k est un entier compris entre 50 et 5000 et f 
10 est exprim6 en Hz ; 

- I'onde ultrasonore focalisee emise au cours de 
I'etape d' excitation (a) presente une frequence comprise 
entre 0^.5 et 15 MHz et est emise pendant une succession de 
periodes d*" emission separees par des periodes de repos, les 

15 periodes d' emission se succedant a une cadence comprise en- 
tre 10 et 1000 emissions par seconde ; 

- I'onde ultrasonore focalisee emise au cours de 
I'^tape d'' excitation (a) est une coiabinaison lineaire (no- 
tamment une somme) de deux signaux monochromatiques de fre- 

20 quences respectives fl et f2 telles que 20 Hz ^ | f 1 - f 2 | ^ 
1000 Hz ; 

I'onde ultrasonore focalisee 6mise au cours de 
l'6tape d' excitation (a) est focalisee siraultan^ment en 
plusieurs points ; 

25 - l'6tape (c) de traitement d' images est suivie 

(immediatement ou non) par une 6tape (d) de cartographie au 
cours de laquelle, a partir d' une evolution du parametre de 
mouvement au cours du temps, on calcule au moins un parame- 
tre de propagation de I'onde de cisaillement en au moins 

30 certains points du champ d' observation pour determiner ain- 
si une cartographie dudit parametre de propagation dans le 
champ d' observation ; 

le parametre de propagation de I'onde de cisail- 
lement qui est calcule au cours de 1 ' etape de cartographie 
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est choisi parmi la vitesse des ondes de cisaillement^ le 
module de cisaillement, le module d' Young, 1' attenuation 
des ondes de cisaillement , I'elasticite de cisaillement , la 
viscosite de cisaillement et le temps de relaxation mecani- 
5 que ; 

on repete successivement les etapes (a) a (d) en 
emettant des ondes ultrasonores differentes au cours des 
etapes d' excitation (a) successives, puis on combine les 
differentes cartographies obtenues au cours des etapes de 

10 cartographie (d) successives, pour calculer une cartogra- 
phie synthetique du champ d' observation. 

Par ailleursr 1' invention a egalement pour objet un 
dispositif d'imagerie utilisant des ondes de cisaillement 
pour observer un milieu visco^lastique diffusant qui 

15 contient des particules ref lechissant les ondes ultrasono- 
res de compression, ce dispositif comprenant un reseau de 
transducteurs commandes ind^pendairanent les uns des autres 
par au moins une unite centrale 61ectronique adaptee pour : 

- faire generer au moins une onde 61astique de ci- 
20 saillement dans le milieu viscoelastique, 

observer la propagation de I'onde de cisaille- 
ment simultan6ment en une multitude.. de--i?ointc d'uu- champ 
d' observation dans le milieu viscoelastique^ en faisant 
emettre dans le milieu viscoelastique, par ledit reseau de 

25 transducteurs, une succession de tirs d' ondes ultrasonores 
de compression non focalisees, a une cadence d'au moins 500 
tirs par seconde, puis en faisant detecter en temps reel 
par ledit reseau de transducteurs, et pour enregistrer en 
temps reel, des signaux acoustiques regus du milieu viscoe- 

30 lastique, comprenant les echos gen^res par les ondes ultra- 
sonores de compression non focalisees en interagissant avec 

les particu l£L^_C€^ f . l 6 ohj rf^.'S-a- n t, b^M rm rlit- rm l t g^n Tri or-r^^l j g^ir ^ 

- et traiter les signaux acoustiques successifs 
re^us du milieu viscoelastique pour determiner des images 
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de propagation successives de l*onde de cisail lament, puis 
determiner au raoins un parametre de mouvement du milieu 
visco61astique en dif f erents points du champ d' observation, 
caracteris6 en ce que l'unit6 centrale electronique est 
5 adaptee pour faire gen6rer I'onde 61astique de cisaillement 
en faisant emettre au moins une onde ultrasonore focalisee 
dans le milieu viscoelastique par ledit reseau de transduc- 
teurs, la focalisation et la chronologie de ladite onde ul- 
trasonore focalisee, ainsi que la chronologie desdites on-- 
10 des ultrasonores non focalis6eSr etant adapt ees pour que 
lesdites ondes ultrasonores non focalisees parviennent dans 
le champ d' observation lors de la propagation de I'onde de 
cisaillement dans ce champ d' observation . 

D'autres caract^ristiques et avantages de 1 ' inven- 
15 tion apparaitront au cours de la description suivante d^ une 
de ses formes de realisation, donnee a titre d'exemple non 
limitatif, en regard du dessin joint. 

Sur le dessin, la figure 1 est une vue schematique 
d'un dispositif d'imagerie par ondes de cisaillement selon 
20 une foarrae de realisation de 1' invention. 

Le dispositif d'imagerie 1 represents sur la figure 
1 est destine a etudier la propagation des ondes elastiques 
de cisaillement dans un milieu viscoelastique 2 qui est 
diffusant vis k vis des ondes ultrasonores de compression, 
25 et qui peut §tre par exemple : 

~ un corps inerte, notamment dans le cas du 
controle de quality pour des applications industrielles, 

ou un corps vivant, par exemple une partie du 
corps d'un patient, dans le cas des applications medicales. 
30 Ces mouvements sont suivis par exemple au moyen 

d'un micro-ordinateur 4 (comprenant au moins une interface 
d' entree 4a telle qu^un clavier ou autre, et une interface 
de sortie telle qu'un ecran ou autre) ou toute autre unit6 
centrale 61ectronique, qui fait envoyer dans le milieu 2, a 
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partir de sa surface exterieure 3^ des ondes ultrasonores 
de compression qui interagissent avec les particules diffu- 
santes 5 contenues dans le milieu 1^ lesquelles particules 
sont ref lechissantes pour les ondes ultrasonores de com- 
5 pression. Les particules 5 peuvent etre constituees par 
toute heterog^neite du milieu 1, et notamment, lorsqu'il 
s'agit d'une application medicale, par des particules de 
collag^ne presentes dans les tissus humains (ces particules 
forment sur les images echographiques des points connus 

10 sous le terme "speckle"). 

Pour observer la propagation de I'onde de cisaille- 
ment, on utilise une sonde ultrasonore 6 disposee contre la 
surface ext6rieure 3 du milieu observe 1. Cette sonde en- 
voie, selon un axe des impulsions d' ondes ultrasonores 

15 de compression du type de celles couramment utilisdes en 
echographie, ^ une frequence comprise par exemple entre 0,5 
et 100 MHz et de preference entre 0,5 et 15 MHz, par exem- 
ple de 1 ' ordre de 4 MHz . 

La sonde ultrasonore 6 est constituee par un r6seau 

20 de n transducteurs ultrasonores Tl, T2, Ti, Tn, n 

6tant un nombre entier n au moins egal ^ 1. 

Cette sonde 6 peut se presenter p^.r- e^crr.ple sou's la" 
forme d'une barrette lineaire pouvant comprendre par exem- 
ple n = 128 transducteurs align^s selon un axe Y perpendi- 

25 culaire k I'axe X. Mais la sonde en question pourrait etre 
egalement un reseau bidimensionnel (plan ou non) de trans- 
ducteurs . 

Les transducteurs Tl, T2, ... Tn sont commandes in- 
dependamment les uns des autres par le micro-ordinateur 4, 
30 eventuellement par 1 ' intermediaire d'une unite centrale CPU 
qui est contenue par exemple dans une bale electronique 7 

reliee par u n r.AhJ,^— sa^^^-e— w- 3h ^^-onde h*: Les transducteurs 

Tl~Tn peuvent ainsi emettre selectivement : 

- soit une onde ultrasonore de compression "plane" 
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(c'est a dire en 1' occurrence une onde dont le front d'onde 
est rectiligne dans le plan X, Y) ou tout autre type d'onde 
non focalisee eclairant 1' ensemble du champ d' observation 
dans le milieu 2^ par exemple une onde g^neree en faisant 
5 emettre des signaux acoustiques aleatoires par les diff6- 
rents transducteurs Tl~Tn, 

- soit une onde ultrasonore de compression focali- 
see en un ou plusieurs points du milieu 2. 

Pour observer la propagation de I'onde de cisaille- 
10 ment dans le milieu 2, on procede en plusieurs etapes suc- 
cessives : 

(a) une etape d' excitation au cours de laquelle le 
micro-ordinateur 4 fait generer une onde elastique de ci- 
saillement dans le milieu viscoelastique 2^ en faisant 

15 emettre au moins une onde ultrasonore focalisee dans- le mi- 
lieu viscoelastique par la sonde 6, 

(b) une etape d' observation au cours de laquelle on 
observe la propagation de I'onde de cisaillement simultane- 
ment en une multitude de points du champ d' observation dans 

20 le milieu viscoelastique 2^ cette etape d' observation com- 
prenant les sous-6tapes suivantes : 

(bl) le micro-ordinateur 4 fait emettre dans 
le milieu viscoelastique, par la sonde 6 une succession de 
tirs d'ondes ultrasonores de compression non focalis6es a 

25 une cadence d'^au moins 500 tirs par seconde (la focalisa- 
tion et la chronologie de l^onde ultrasonore focalisee 
emise k I'^tape (a), ainsi que la chronologie desdites on- 
des ultrasonores non focalisees, sont adaptees pour qu'au 
moins certaines desdites ondes ultrasonores non focalis6es 

30 parviennent dans le champ d' observation lors de la propaga- 
tion de I'onde de cisaillement dans ce champ d' observation, 
pour au moins certaines emissions d'onde ultrasonore non 
focalisee) , 

(b2) le micro-ordinateur 4 fait detecter par la 
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sonde 6 et enregistrer en temps reel des signaux acousti- 
ques reg:us du milieu viscoelast ique 2^ comprenant les echos 
generes par les ondes ultrasonores de compression non foca- 
lisees en interagissant avec les particules ref lechissantes 
5 5 du milieu viscoelastique, ces 6chos correspondant (direc- 
tement ou indirectement ) a des images successives de depla- 
cement du milieu viscoelastique, 

(c) et au moins une etape de traitement au cours de 
laquelle : 

10 (cl) le micro-ordinateur 4 traite les signaux 

acoustiques successifs regus du milieu viscoelastique 2 au 
cours de la sous-etape (b2) pour determiner des images de 
propagation successives, 

(c2) et le micro-ordinateur 4 determine au moins 

15 un parametre de mouvement du milieu viscoelastique 2 en 
diff^rents points du champ d' observation. 

L'onde ultrasonore focalis^e emise au cours de 
1' etape d' excitation (a) peut Stre une onde monochromatique 
de frequence f comprise entre 0,5 et 15 MHz, par exemple 

20 environ 4 MHz, emise pendant une duree k/f en secondes, oti 
k est un entier compris entre 50 et 5000 (par exemple de 
I'ordre de 500) et f est exprime en -Evsntusilerutiiic, une 

telle onde peut etre emise pendant une succession de perio- 
des d' emission s^parees par des p^riodes de repos, les pe- 

25 riodes d' emission se succedant a une cadence comprise entre 
10 et 1000 emissions par seconde. 

En variante, I'onde ultrasonore focalisee emise au 
cours de 1' etape d' excitation (a) est une combinaison li- 
neaire (notamment une somme) de deux signaux monochromati- 

30 ques de frequences respectives fl et f2 telles que 20 Hz <^ 
I f 1 ~ f 2 I ^ 1000 Hz, ce qui produit une onde modulee en am- 

Plitude a ver nn e^-^^eq^ iiL ^ H d^ - 'i ffoduJ at -i nn | f 1 - f 2 [ . 

Par ailleurs, 1' onde ultrasonore focalisee emise au 
cours de 1' etape d' excitation (a) peut eventuellement etre 
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focalisee, simultan^ment* ou non, en plusieurs points, de 
fa?on que I'onde de cisaillement gen6ree pr§sente la forme 
souhait^e (par exemple, on peut ainsi gen§rer une onde de 
cisaillement plane, ou au contraire une onde de cisaille- 
5 ment focalisee) et illumine les zones souhait^es du milieu 
2. 

Au cours de I'etape (bl) , qui peut durer par exem- 
ple moins d'une seconde, on peut 6mettre les ondes ultraso- 
nores de compression non focalisees, a une cadence comprise 

10 entre 500 et 10 000 tirs par seconde et de preference com-- 
prise entre 1000 et 5 000 tirs par seconde (cette cadence 
est limitee par le temps d ' aller-retour de I'onde de com- 
pression dans le milieu 2, done par I'epaisseur du milieu 2 
dans la direction X : il faut en effet que tous les echos 

15 generes par I'onde de compression aient ete regus par la 
sonde 6 avant d'envoyer une nouvelle onde de compression). 

Chaque onde ultrasonore de compression non focali- 
see se propage dans le milieu 2 avec une vitesse de propa- 
gation beaucoup plus ^levee que les ondes de cisaillement 

20 (par exemple de I'ordre de 1500 m/s dans le corps humain) , 
et interagit avec les particules r6f 16chissantes 5, ce qui 
genere des 6chos ou autres perturbations analogues du si- 
gnal, connus en soi sous le nom de "bruits de speckle" dans 
le domaine de 1 ' echographie. 

25 Ces "bruits de speckle" sont captes par les trans- 

duct eurs Tl, ... ,Tn au cours de la sous-etape (b2) , apres 
chaque tir d'onde ultrasonore de compression non focalisee. 
Le signal sij (t) ainsi capt6 par chaque transducteur Ti 
aprds le tir n*" j est tout d'abord 6chantillonn6 a haute 

30 frequence (par exemple de 30 k 100 MHz) et numerise en 
temps reel (par exemple sur 12 bits) par un echantillonneur 
appartenant a la bale 7 et reli^ ci ce transducteur, respec- 
tivement El, E2, ... En- 

Le signal sij (t) ainsi echantillonne et numerise 
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est ensuite memorise, 6galement en temps r6el,. dans une me- 
moire Mi appartenant k la baie 7 et propre au transducteur 
Ti. 

Chaque memoire Mi presents par exemple une capacite 
5 de I'ordre de 128 Mo^ et contient 1 ' ensemble des signauK 
sij (t) regus successivement pour les tirs j = 1 ^ p. 

En temps differe, apres la memorisation de tous les 
signaux sij (t) correspondant a une meme propagation d'onde 
de cisaillement, 1^ unite centrale CPU fait retraiter ces 
10 signaux par un circuit sommateur S appartenant a la baie 7 
(ou bien elle effectue elle-mSme ce trait ement, ou encore 
ledit traitement peut etre effectu6 dans le micro- 
ordinateur 4), par un processus classique de formation de 
voies correspondant S la sous-etape (cl) • 
15 On genere ainsi des signaux Sj (x, y> qui corres- 

pondent chacun a 1' image du champ d' observation apres le 
tir j (dans le cas ou I'onde ultrasonore non focalisee 
est une onde plane) . 

Par exemple, on peut determiner un signal Sj (t) 
20 par la formule suivante : 

Sj\t)^^a,ix,y).sim^^ -~ - 

- - - • - - * 1=1 

o\a : 

sij est le signal brut pergu par le transducteur 
n*" i apres le tir d' onde ultrasonore de compression n** j, 
25 - t(x^y) est le temps mis par I'onde ultrasonore 

de compression pour atteindre le point du champ 
d' observation de coordonn6es (x^y), avec t = 0 au debut du 
tir n** j, 

di{x,y) est la distance entre le point du champ 
30 d' observati on de cQQrdo nnp>g>c^ (^rV) i^-^fcL^.i=>->:>A^^-^^ — rr^'"T7 — 
ou une approximation de cette distance^ 

V est la Vitesse moyenne de propagation des on- 
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des acoustiques ultrasonores de compression dans le milieu 

viscoelastique observe . 

et ai(x,y) est un coefficient de ponderation 

tenant compte de lois d' apodisation (en pratique, on pourra 
5 dans de nombreux cas considerer que ai(x^y)=l). 

La formule ci-dessus s' applique mutatis mutandis 

lorsque le champ d' observation est a 3 dimensions (reseau 

bidimensionnel de transducteurs) , en rempla<jant les coor- 

donnees spatiales (x,y) par (x,y,z). 
10 Apres I'etape eventuelle de formation de voies, 

1' unite centrale CPU memorise dans una memoire centrale M 

appartenant a la bale 7 les signaux d'' images Sj(x,y) ou 

Sj (x) ou Sj(x,y,z), qui correspondent chacun au tir n** " j . 

Ces signaux peuvent egalement etre memorises dans le micro- 
15 ordinateur 4 lorsqu'il effectue lui-m§me le traitement 

d' image . 

Ces images sont ensuite traitees en temps differ§ ^ 
la sous-etape (c2) , par correlation et avantageusement par 
intercorrelation soit deux it deux, soit de preference avec 

20 une image de reference qui peut 6tre : 

- soit une image de d^placement precedeironent de- 
termines comme expliqu6 ci-dessus et utilisee comme image 
de reference pour les images de deplacement ulterieures (ou 
pour un nombre limite d' images de d§placement ulterieures, 

25 par exemple 30 images de deplacement) , 

soit determin^e au cours d^une etape prelimi- 
naire d' observation initiale (aO) , comme les images succes- 
sives de deplacement susmentionnees, en faisant emettre une 
ou plusieurs ondes ultrasonores non focalisees par la sonde 

30 6 avant l'6tape d' excitation (a) qui g6nere I'onde de ci- 
saillement (lorsque plusieurs ondes ultrasonores de com- 
pression non focalisees sont ainsi emises avant la phase 
d' excitation, on enregistre des echos generes par chaque 
onde ultrasonore de compression non focalisee en interagis- 
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sant avec les particules ref lechissantes du milieu viscoe- 
lastique, ces echos correspondant a plusieurs images preli- 
minaires successives du milieu viscoelastique, et on deter- 
mine ladite image initiale du milieu viscoelastique en com- 
5 binant lesdites images preliminaires successives et notam-- 
ment en moyennant les valeurs des pixels desdites images 
preliminaires) • 

L ' intercorr61ation susmentionnee peut §tre realises 
par exemple dans un circuit §lectronique specialise DSP ap- 
10 partenant k la bale 7, ou etre programmee dans 1' unite cen- 
trals CPU ou dans le rtiicro-ordinateur 4. 

Au cours de ce processus d' intercorrelation, on 
maximise une fonction d' intercor relation 

<Sj (x,y) ,Sj + l (x,y)> afin de determiner le d6placement subi 
15 par chaque particule 5 donnant lieu a un echo ultrasonore. 

Des exemples de tels calculs d' intercorrelation 
sont donnas dans l'6tat de la technique, notarament par 
O'Donnell et al. ("Internal displacement and strain imaging 
using speckle tracking"/ IEEE transactions on ultrasonic, 
20 f erroelectrics, and frequency control, vol. 41, n*" 3, mai 
1994, p, 314-325) et par Ophir et al. ( "Elastography : a 
quantitative method for imagincr the .el'?^5?ti^it-y -of • bio'loyi- 
cal tissues". Ultrasonic imag., vol* 13, p-111-134, 1991). 

On obtient ainsi un ensemble de vecteurs deplace- 
25 ments engendres par les ondes de cisaillement en cha- 

que position ^ du milieu 2 sous I'effet de I'onde de ci- 
saillement (ces vecteurs deplacements peuvent eventuelle- 
ment etre reduits a une seule composante dans 1' exemple 
considere ici) . 

30 Get ensemble de vecteurs deplacements est stocke 

dans la memoire M ou dans le micro-ordinateur 4 et peut par 

exemple pfrp vi s4^-3.4-a^ — nu L dimue 'nt au moyen de 1^ ecran 4a du 

micro-ordinateur, sous la forme d'un film ralenti oix la va- 
leur des deplacements est illustree par un parametre opti- 
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que tel que par un niveau de gris ou par un niveau chroma- 
tique . 

On visualise ainsi parf aitement les differences de 
propagation de I'onde de cisaillement entre les zones de 
5 caract^ristiques differentes du milieu 2^ par exemple les 
tissus sains et les tissus cancereux dans le cas d'une 
application mfedicale. 

Ce film de propagation de I'onde de cisaillement 
est en outre superposable avec une image echographique 

10 classique^ qui peut etre generee egalement par le disposi-- 
tif 1 decrit ci-dessus. 

Par ailleurs; il est egalement possible de calculer 
non pas les d^placements de chaque point du milieu observe 
2f mais les deformations du milieu 2, c'est a dire des vec- 

15 teurs dont les composantes sont les derivees des composan- 
tes des vecteurs deplacements respectivement par rapport 
aux variables d'espace (coordonnees selon X, Y dans 1' exem- 
ple considere) • Ces vecteurs de deformation sont utilisa- 
bles comme les vecteurs deplacements pour visualiser clai- 

20 rement la propagation de I'onde de cisaillement sous la 
forme d'un film, et presentent en outre l^avantage de s'af- 
franchir des deplacements de la sonde 6 par rapport au mi- 
lieu observe 2 . 

A partir des champs de deplacements ou de deforma- 

25 tions, le micro-ordinateur 4 peut avantageusement proc^der 
ensuite k une 6tape de cartographic (d) au cours de la- 
quelle, a partir de 1* Evolution du parametre de mouvement 
(deplacement ou deformation) au cours du temps dans le 
champ d' observation X, Y (ou X, Y, Z dans le cas d'un re- 

30 seau bidimensionnel de transducteurs) , on calcule au moins 
un parametre de propagation de I'onde de cisaillement, soit 
en certains points du champ d' observation choisis par 
1 ' utilisateur a partir du micro-ordinateur 4, soit dans 
tout le champ d' observation. 



25 
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Le parametre de .propagation de I'onde de cisaille- 
ment qui est calculi au cours de I'etape de cartographie 
est choisi par exemple parmi : la vitesse Cs des ondes de 
cisaillement, le module de cisaillement [i on le module 
5 d' Young E=3|li^ 1 ' attenuation a des ondes de cisaillement, 
I'elasticite de cisaillement jil, la viscosite de cisaille- 
ment |a2, ou le temps de relaxation mecanique ts des tissus. 

Par exemple, on peut calculer en differents points 
du champ d' observation : 
10 - la valeur de la celerite Cs de I'^onde de cisail- 

lement, qui donne acces a la durete des tissus, 

la valeur du temps de relaxation m6canique ts 
des tissus, caracteristique de la viscosite locale du mi- 
lieu, 

15 Pour cela, on utilise 1' Equation de propagation (1) 

suivante, a laquelle obeissent les d6placements e^- 
gendr^s par les ondes de cisaillement en chaque position ^ 
du milieu : 

P^5^ = c/(l-.r,|-).V^«(F.O 
20 (1) 

oti p est la density des tissue, xS est. .le temp?, de rsiaxa- 

tion m6canique des tissus et cS est la cel6rit6 de I'onde 

de cisaillement, directement relive au module d' Young E des 

tissus par la relation : 



= 1^ 

VP 



(2) 



Dans le domaine de Fourier, 1' Equation d'onde (1) 
ci-dessus peut s'ecrire : 



30 
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est la transforme de Fourier du champ de deplace- 
ment ^^^^^ mesure en chaque point et ^^(^»^)est la trans- 
f orniee de Fourier du laplacien spatial de ce champ "(^'^) , 
Etant donn6 que coxS « 1, on a done acc6s h une expression 
simplifiee : 

Cs ^co p 



AU(r,a>) (4 J 

(D AU(r,a)) 



ou est la phase de la variable complexe x, 

Les fonctions ^('^'^)et ^U(r,a>) ^tant connues en chaque 
point de 1' image echographique^ il est done possible de 
mesurer en tout point de I'espace le module d' Young et le 

15 temps de relaxation m§canique des tissus, en 6tablissant 
ainsi une cartographie de ces deux parametres • 

De plus^ les equations (4) et (5) etant vraies a 
chaque frequence le calcul de cS et xS peut etre avanta- 
geusement moyenne sur 1' ensemble de la bande des frequences 

20 portees par I'onde de cisaillement ^ ameliorant ainsi forte- 
ment la qualite de la cartographie realisee. A cet effete 
on peut utiliser les formules suivantes : 



= w^p — ^ - da> 



(6) 



(7) 



ou coO et G) 1 sont les frequences minimales et maximales que 
30 porte I'onde de cisaillement 
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Le mode de calcul serait le mSme en utilisant non 
plus les d^placements, mais les deformations du milieu ob- 
serve 2 . 

Par ailleurs, on peut avantageusement etablir suc^ 
5 cessivement plusieurs cartographies des parametres de pro-- 
pagation souhaites, par exemple cS et xS, en generant sue- 
cessivement des ondes de cisaillement dif f erentes;. obtenues 
par exemple en emettant des ondes ultrasonores de compres- 
sion focalisees successivement en plusieurs points ou ayant 
10 des formes d' ondes differentes. On peut ensuite combiner 
les differentes cartographies obtenues^ par exemple en les 
moyennant, de fagon a obtenir une cartographie synthetique 
plus riche et plus precise. 
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REVENDXCATIONS 



1. Proc6de d'imagerie utilisant des ondes de ci- 
saillement pour observer un milieu visco61astique diffusant 
(2) qui contient des particules (5) ref 16chissant les ondes 
ultrasonores de compression, proc§d6 comprenant : 

(a) une 6tape d' excitation au cours de laquelle on 
genere une onde elastique de cisaillement dans le milieu 
viscoelastique (2), 

(b) une etape d' observation au cours de laquelle on 
observe la propagation de I'onde de cisaillement simultan^- 
ment en une multitude de. points d'un champ d' observation 
dans le milieu, viscoelastique (2), cette etape d'observa- 
tion comprenant les sous-etapes suivantes : 

(bl) faire emettre dans le milieu viscoelas- 
tique (2) r par un reseau de transducteurs (6) commandes in- 
dependamment les uns des autres^ une succession de tirs 
d' ondes ultrasonores de compression non focalisees a une 
cadence d' au moins 500 tirs par seconde, 

(b2) faire detecter et enregistrer en temps reel 
des signaux acoustiques regus du milieu viscoelastique (2), 
comprenant les echos generes par les ondes ultrasonoa;'es de 
compression non focalisees en interagissant avec les parti- 
cules ref l§chissantes (5) dudit milieu viscoelastique, 

(c) et au moins une etape de traitement au cours de 
laquelle : 

(cl) on traite les signaux acoustiques succes- 
sifs regus du milieu viscoelastique (2) au cours de la 
sous-etape (b2) pour determiner des images de propagation 
successives de I'onde de cisaillement, 

(c2) et on determine au moins un param^tre de 
mouvement du milieu viscoelastique (2) en diff^rents points 
du champ d' observation, 

caracterise en ce qii' au cours de I'^tape d' excitation (a). 
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on fait generer I'onde 61astique de cisaillement en faisant 
emettre au moins une onde ultrasonore focalisee dans le mi- 
lieu visco§lastique (2) par ledit r§seau de transducteurs 
(6), la focalisation et la chronologie de ladite onde ul- 
trasonore focalisee, ainsi que la chronologie desdites on- 
des ultrasonores non focalis6es, etant adapt6es pour qu ' au 
moins certaines desdites ondes ultrasonores non focalisees 
parviennent dans le champ d' observation lors de la propaga- 
tion de I'onde de cisaillement dans ce champ d' observation, 
pour au moins certaines Emissions d'onde ultrasonore non 
focalis6e. 

2. Precede selon la revendication 1, dans lequel, 
au cours de la sous-etape (b2) , pour determiner ledit para- 
metre de mouvement, on compare plusieurs images successives 
de propagation avec une m§me image de r6f6rence du milieu 
viscoelastique (2), cette image de r6f6rence etait d^termi- 
n6e en tirant au moins une onde ultrasonore de compression 
non focalisee dans ledit milieu viscoelastique puis en d6- 
tectant et en enregistrant des 6chos g^nires par ladite 
onde ultrasonore de compression non focalisee lorsqu'elle 
interagit avec les particules r6f 16chissantes (5) du milieu 



3. Precede selon la revendication 2, dans lequel 
l'6tape (a) est pr6cedee d'une etape d' observation initiale 
(aO) au cours de laquelle on tire au moins une onde ultra- 
sonore de compression non focalis6e puis on detecte et on 
enregistre des echos g6n6res par ladite onde ultrasonore de 
compression non focalisee en interagissant avec les parti- 
cules reflechissantes (5) du milieu viscoelastique, ces 
6chos correspondant a une image initiale du milieu vis^coe- 
lastique, et au cours de la sous-etape (b2) ladite image 
initiale constitue ladite image de reference pour traiter 
au moins certaines des images successives de d6placement. 

4. Precede selon la revendication 3, dans lequel. 
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au cours de I'^tape d' observation initiale (aO) , on tire 
successivement plusieurs ondes ultrasonores de compression 
non focalisees puis on detecte et on enregistre des echos 
generes par chaque onde ultrasonore de compression non fo- 
5 calisee en interagissant avec les particules ref lechissan- 
tes (5) du milieu viscoelastique, ces echos correspondant a 
plusieurs images successives du milieu viscoelastique^ et 
on determine ladite image initiale du milieu viscoelastique 
en combinant lesdites images successives, 
10 5. Precede selon I'une quelconque des revendica- 

tions pr6cedentes, dans lequel ledit parametre de mouvement 
est un deplacement du milieu viscoelastique (2) . 

6. Proc6d6 selon I'une queldonque des revendica- 
tions precedentes, dans lequel I'onde ultrasonore focalisee 

15 emise au cours de I'etape d' excitation (a) presents une 
frequence f comprise entre 0,5 et 15 MHz et est 6mise pen- 
dant une duree k/f en secondes, oil k est un entier compris 
entre 50 et 5000 et f est exprime en Hz. 

7. Precede selon I'une quelconque des revendica- 
20 tions 1 a 5, dans lequel l^'onde ultrasonore focalis6e emise 

au cours de I'etape d' excitation (a) pr^sente une frequence 
comprise entre 0,5 et 15 MHz et est emise pendant une suc- 
cession de p6riodes d' emission s6parees par des p6riodes de 
repos, les pferiodes d' Emission se succ^dant ci une cadence 
25 comprise entre 10 et 1000 emissions par seconde. 

8. Precede selon I'une quelconque des revendica- 
tions 1 a 5, dans lequel I'onde ultrasonore focalisee emise 
au cours de l''etape d^ excitation (a) est une combinaison 
lineaire (notamment une somme) de deux signaux monochroma- 

30 tiques de frequences respectives fl et f2 telles que 20 Hz 
< |fl - f2 I <> 1000 Hz. 

9. Precede selon I'une quelconque des revendi ca- 
tions precedentes, dans lequel I'onde ultrasonore focalisee 
emise au cours de I'etape d' excitation (a) est focalisee 
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sixnultanement en plusieurs points. 

10. Precede selon I'une quelconque des revendica- 
tions precedentes^ dans lequel 1 ' etape (c) de traitement 
d' images est suivie par une etape (d) de cartographie au 
5 cours de laquelle, a partir d' une Evolution du parametre de 
mouvement au cours du temps, on calcule au moins un parame- 
tre de propagation de I'onde de cisaillement en au moins 
certains points du champ d' observation pour determiner ain- 
si une cartographie dudit parametre de propagation dans le 

10 champ d' observation • 

11- Precede selon I'une quelconque des revendica- 
tions precedentes, dans lequel le parametre de propagation 
de I'onde de cisaillement qui est calcule au cours de 
l'6tape de cartographie (d) est choisi parmi la vitesse des 

15 ondes de cisaillement, le module de cisaillement, le module 
d' Young, 1' attenuation des ondes de cisaillement, 
I'elasticite de cisaillement, la viscosite de cisaillement 
et le temps de relaxation m6canique. 

12. Precede selon la revendication 11, dans lequel 
20 on repete success ivement les etapes (a) a (d) en emettant 

des ondes ultrasonores focalisees differentes au cours des 
etapes d' excitation (a) successives, puis on combine les 
differentes cartographies obtenues au cours des etapes de 
cartographie (d) successives, pour calculer une cartogra- 
25 phie synthetique du champ d' observation . 

13. Dispositif d'imagerie pour la mise en ceuvre 
d'un procede selon I'une quelconque des revendications pre- 
cedentes, utilisant des ondes de cisaillement pour observer 
un milieu viscoelastique diffusant (2) qui contient des 

30 particules (5) ref lechissant les ondes ultrasonores de com- 
pression, ce dispositif comprenant un reseau de transduc- 

^Pi^^T- ' *^ — C£4 — r . ommande s — independamment — i^s — un s d e s — autrcp par 

au moins une unite cent rale eiectronique (4, CPU) adaptee 
pour : 
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faire g6n6rer au moins une onde ^lastique de ci- 
saillement dans le milieu viscoelastique (2), 

- observer la propagation de I'onde de cisaille- 
ment simultanement en une multitude de points d'un champ 
5 d' observation dans le milieu viscoelastique (2), en faisant 
emettre dans le milieu viscoelastique, par ledit reseau de 
transducteurs (6), une succession de tirs d' ondes ultraso- 
nores de compression non focalisees, ^ une cadence d' au 
moins 500 tirs par seconde, puis en faisant detecter en 

10 temps reel par ledit reseau de transducteurs (6), et pour 
enregistrer en temps reel, des signaux acoustiques regus du 
milieu viscoelastique (2), comprenant les echos gen^res par 
les ondes ultrasonores de compression non focalisees en in- 
teragissant avec les particules ref lechissantes (5) dudit 

15 milieu viscoelastique, 

et traiter les signaux acoustiques successifs 
iregus du milieu viscoelastique (2) pour determiner des ima- 
ges de propagation successives de I'onde de cisaillement , 
puis determiner au moins un param^tre de mouvement du mi- 

20 lieu viscoelastique (2) en differents points du ch^mp 
d' observation, 

caracterlse en ce que 1' unite cent rale electronique (4) est 
adaptee pour faire g6nerer I'onde elastique de cisaillement 
en faisant 6mettre au moins une onde ultrasonore focalisee 

25 dans le milieu viscoelastique par ledit reseau de transduc- 
teurs (6), la focalisation et la chronologie de ladite onde 
ultrasonore focalisee, ainsi que la chronologie desdites 
ondes ultrasonores non focalisees, etant adaptees pour que 
lesdites ondes ultrasonores non focalisees parviennent dans 

30 le champ d' observation lors de la propagation de I'onde de 
cisaillement dans ce champ d' observation. 
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